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Abstract 

We discuss the production of massive relic coherent gravitons in a 
particular class of f{R) gravity which arises from string theory and their 
possible imprint in Cosmic Microwave Background. In fact, in the very 
early universe these relic gravitons could have acted as slow gravity waves. 
They may have then acted to focus the geodesies of radiation and matter. 
Therefore, their imprint on the later evolution of the universe could ap- 
pear as filaments and domain wall in the Universe today. In that case, the 
efi^ect on Cosmic Microwave Background should be analogous to the effect 
of water waves, which, in focusing light, create optical caustics which are 
commonly seen on the bottom of swimming pools. We analyze this im- 
portant issue by showing how relic massive GWs perturb the trajectories 
of Cosmic Microwave Background photons (gravitational lensing by relic 
GWs). 

The consequence of the type of physics discussed is outhned by illus- 
trating an amplification of what might be called optical chaos. 



1 Introduction 

Modified gravity currently obtains a lot of attention from the scientific commu- 
nity. The main reason is the remarkable issue that it enables a description of 
early-time inflation as well as late-time acceleration epoch (Dark Energy) in a 
unified way. 

In recent years, superstring/M theory caused a lot of interest about higher 
order gravity in more than 4 dimensions [T] . These models work in the effective 
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low-energy action of superstring theory pQ [5] . Within the classical framework, 
they have to be inserted among the class of the so-called f{R) theories of gravity 
(for a recent review see [3]). 

Motivations for a potential extension of Einstein's general relativity (GR) [J] 
are various. First of all, as distinct from other field theories, like the electromag- 
netic theory, GR is very difficult to quantize. This fact rules out the possibility of 
treating gravitation like other quantum theories and precludes the unification of 
gravity with other interactions. At the present time, it is not possible to realize 
a consistent quantum gravity theory which leads to the unification of gravita- 
tion with the other forces. One of the most important goals of modern physics 
is to obtain an unified theory which could, in principle, show the fundamental 
interactions as different forms of the same symmetry. Considering this point of 
view, today one observes and tests the results of one or more breaks of symme- 
try. In this way, it is possible to say that we live in an unsymmetrical world [S] . 
In the last 60 years, the dominant idea has been that a fundamental description 
of physical interactions arises from quantum field theory [5] . In this approach, 
different states of a physical system are represented by vectors in a Hilbert space 
defined in a space-time, while physical fields are represented by operators (i.e. 
linear transformations) on such a Hilbert space. The greatest problem is that 
this quantum mechanical framework is not consistent with gravitation, because 
this particular field, i.e. the metric 5^^, describes both the dynamical aspects 
of gravity and the space-time background [S]. In other words, one says that 
the quantization of dynamical degrees of freedom of the gravitational field is 
meant to give a quantum-mechanical description of the space-time. This is an 
unequaled problem in the context of quantum field theories, because the other 
theories are founded on a fixed space-time background, which is treated like a 
classical continuum. Thus, at the present time, an absolute quantum gravity 
theory, which implies a total unification of various interactions, has not been ob- 
tained [5] . In addition, GR assumes a classical description of the matter which 
is totally inappropriate at subatomic scales, which are the scales of the early 
Universe [31 [5]. 

In the general context of cosmological evidence, there are also other con- 
siderations which suggest an extension of GR [31 [7]. As a matter of fact, the 
accelerated expansion of the Universe, which is observed today, implies that 
cosmological dynamics is dominated by the so called Dark Energy, which gives 
a large negative pressure. This is the standard picture, in which this new in- 
gredient should be some form of un-clustered, non-zero vacuum energy which, 
together with the clustered Dark Matter, drives the global dynamics. This is 
the so called ^^concordance modeV^ (ACDM) which gives, in agreement with 
the Cosmic Microwave Background Radiation, Large Scale Structure and Su- 
pernovae la data, a good picture of the observed Universe today, but presents 
several shortcomings such as the well known ^' coincidence^^ and ^^Cosmological 
Constant" problems [5]. 

An alternative approach is seeing if the observed cosmic dynamics can be 
achieved through an extension of GR |3l [J. In this different context, it is not 
required to find candidates for Dark Energy and Dark Matter that, untill now. 
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have not been found; only the ^^observed" ingredients, which are curvature and 
baryon matter, have to be taken into account. Then, Dark Energy and Dark 
Matter have to be considered hke pure effects of the presence of an intrinsic 
curvature in the Universe. Considering this point of view, one can think that 
gravity is different at various scales and there is room for alternative theories. 

Note that we are not claiming that GR is wrong. It is well known that, even 
in the context of extended theories of gravity, GR remains the most important 
part of the structure [7] . We are only trying to understand if weak modification 
of could be needed to solve some theoretical and current observational problems. 
In this picture, we also recall that even Einstein tried to modify the framework 
of GR by adding the Cosmological ConstanV^ [3]. In any case. Cosmology 
and Solar System tests show that modifications of GR in the sense of extended 
theories of gravity have to be very weak [3j [7] . 

In principle, the most popular Dark Energy and Dark Matter models can be 
achieved in the framework of extended theories of gravity, i.e. f{R) theories of 
gravity [5] and scalar tensor theories of gravity [7] which are generalizations of 
the Jordan-Fierz-Brans-Dickc Theory [TUl [TTl [T2] . One assumes that geometry 
(for example the Ricci curvature scalar R) interacts with material quantum 
fields generating back-reactions which modify the gravitational action adding 
interaction terms (examples are high-order terms in the Ricci scalar and/or in 
the Ricci tensor and non minimal coupling between matter and gravity). This 
approach enables the modify of the Lagrangian, with respect to the standard 
Einstein- Hilbert gravitational Lagrangian [13], through the addition of high- 
order terms in the curvature invariants (terms like i?^, R°'^Rap, R"^^^^ Rap-^s, 
RDR, RD'^R, in the sense of f{R) Theories [H [7]) and/or terms with scalar 
fields non- minimally coupled to geometry (terms like (j>^R) in the sense of Scalar- 
Tensor Theories 0. 

In the tapestry of /(i?) theories, the higher order terms are physically a 
type of back reaction from geometry acting upon matter which further modifies 
geometry. This is a topological massive gravity which represents a form of 
intrinsic curvature to spacetime. These terms are related to the Bel-Robinson 
tensor [T^ 

vap — R^ ^ aRvapp ^" R^ ^ Rva[ip "^^^ vR ^'^ trRa/S^p- (1) 

Contraction over indices gives the result 

S^g'^^T'^.ap = R'^'^Rpp + R^^^R'^o.pp - (2) 

The physical consequences of this extension to curvature are fairly remarkable. 
The Bel-Robinson tensor is a vacuum curvature VT = 0, and it predicts gravity 
waves (GWs). 
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2 Gravity waves in f{R) theories 



In f{R) gravity the GWs have longitudinal structure (TKTSIIIS], which makes 
them comparable to acoustical waves in a media. The linearized theory of weak 
GWs with a metric perturbation [151 Ull ITT] 



gives a traceless solution in standard GR [T7] 



(3) 



Uh = 0. 



(4) 



The modified gravity results in a terms that acts as a mass, where the wave 
equation is [HI |T2] 



Uh = rn^h. (5) 

The decomposition of the solution gives the standard and /ixx polarization 
modes, the mass introduces a third polarization which is a longitudinal mode 

[Tung. 

Let us consider a string theory setting [I]. The gravitational action is ex- 
panded in powers of a'"i?^" [2], for a the string parameter. The action to 
0(a')is [2] 



5*: 



2k 



L 



(6) 



being L the Lagrangian for everything else compactified on a ZJp-brane. The 
extremization of this gives 



^S* = = / 



2k Sgf^^ ' Sg^^ 



5L 



-gd'^x 



+ 1 



^R + a'R^'^^PR^.^p 



L 



2 s{-g} M„ 



This action derives a modified form of the Einstein field equation 

R 



(7) 



(8) 



Now, let us consider a metric with the form ([3]), which is a classical expres- 
sion. The quadratic term corresponds to quantum corrections on the order of 
the parameter a. We consider this correction as due to fields 0^ so that the 
quantum correction to the metric is 



gtiu ^ Vt^u + hfj^u + 4''^(t)ua, (9) 

where we can regard (ff^cfiya = Sh^^^. These fields are physically a quantum 
correction to the classical gravitational radiation /i^^. In general, these fields 
are quantized fields. In a string theory framework [1] we may define operators 
of the form 



4 



«rn+„«-„, (10) 

m,n— 1 

which is a harmonic oscillator quantization condition compatible with a string 
theory interpretation [T]. The graviton fields are given by the level matching 
condition on a closed string m = 0, which for the first mode n = 1 are states of 
the form 

<^-0,, = a5'y(a;)e'^"a!liy(a;)e''="' (11) 

such that a'la'^_i\0) = constructs the elementary states. 

The connection terms are computed as 

= d,rp€, (12) 

where the field is treated as a vierbein. Now, let us compute the curvature as 

i?^.p = a,<-9p^(;. + K,c^p]^ (13) 

The connection terms arc on the order of the fundamental length a' , which is 
small enough to ignore the second order term. The quantized graviton field may 
then be written in this linearized fashion as 

(14) 

- {dM^Wl, - (9^9,0^)02 + d,cp^^d„4,i - d,<j>';dp<j>2, 

where the last term is zero in a linearized approximation. 

The linearized approximation occurs for long wavelength gravitons. Assume 
the connection term uij^^ = d^cp'^cj)'^ is eigenvalued with a wave number fc" 

= K^^^i (15) 

so the curvature tensor is 

(16) 

= (fcaA:.)0i;0^ - kpKc^i-^cpl = {KK)5h^; - kpkjh^. 
The second order term in the action is then 



R^.uapR'"""' ^ [{kaku)5h^'p - {kpk,)5K]Wk'')5K, - {kPk'')5hp,] = 6k\ (17) 

The term ak'^ is a quartic term in mass, where the string coupling constant 
~ Gn has naturalized units of area. This is an intrinsic curvature in spacetime. 
The string coupling constant is about a' ~ 10^^*^ cm^, which is a small number. 
This also guarantees the viability of the action ([6]) because the theory can pass 
Solar System and Cosmology tests [7]. 
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Is it possible that this mass effect should then become apparent in the labo- 
ratory? The question is, what is the laboratory? The obvious laboratory is the 
Cosmic Microwave Background (CMB). In fact, we recall that relic gravitons 
should have been produced in the Inflationary Era. This is a consequence of 
general assumptions. Essentially it derives from a mixing between basic princi- 
ples of classical theories of gravity and of quantum field theory [THl [IHl HO] ■ The 
strong variations of the gravitational field in the early universe amplify the zero- 
point quantum oscillations and produce relic GWs. It is well known that the 
detection of relic GWs is the only way to learn about the evolution of the very 
early universe, up to the bounds of the Planck epoch and the initial singularity 
[IHl HB HO]- It is very important to stress the unavoidable and fundamental 
character of this mechanism. The model derives from the inflationary scenario 
for the early Universe [T9|, which is tuned in a good way with the WMAP 
data on the CMB (in particular exponential inflation and spectral index w 1) 
pn \n\ . Inflationary models of the early Universe were analyzed in the early 
and middles 1980's [12]. These are cosmological models in which the Universe 
undergoes a brief phase of a very rapid expansion in early times. In this context 
the expansion could be power-law or exponential in time. Inflationary models 
provide solutions to the horizon and flatness problems |19] and contain a mech- 
anism which creates perturbations in all flelds [THl HO] • Important for our goals 
is that this mechanism also provides a distinctive spectrum of relic GWs. The 
GWs perturbations arise from the uncertainty principle and the spectrum of 
relic GWs is generated from the adiabatically-amplified zero-point fluctuations 
[181I20]. 

Relic gravitons can be characterized by a dimensionless spectrum [HI [20] 



where 



. (19) 

is the (actual) critical density energy, pc of the Universe, Hq the actual value 
of the Hubble expansion rate and dpgw the energy density of relic GWs in the 
frequency range / to f + df. 

In standard inflationary model the spectrum is flat over a wide range of 
frequencies, see [HI [20] and figure 1. The more recent value for the flat part of 
the spectrum that arises from the WMAP data can be found in [20] . 

^gUf) < 9 X 10-13 (20) 
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Figure 1: adapted from ref. [5D] 



The spectrum of relic scalar GWs in inflationary models is flat over a wide range 
of frequencies. The horizontal axis is logj^Q of frequency, in Hz. The vertical 
axis is logj^Q ilgsw ■ The inflationary spectrum rises quickly at low frequencies 
(wave which re-entered in the Hubble sphere after the Universe became matter 
dominated) and falls off above the (appropriately redshiftcd) frequency scale 
fmax associated with the fastest characteristic time of the phase transition at 
the end of inflation. The amplitude of the flat region depends only on the energy 
density during the inflationary stage; we have chosen the largest amplitude 
consistent with the WMAP constrains on scalar perturbations. This means 
that at LIGO and LISA frequencies, r2g„,(/)/i^g„ < 9 * 10^^^ 

Based on the weakness of the signal, it will be very difficult to detect relic 
gravitons on Earth, but a potential detection could be, in principle, realized 
with LISA |15j . However, the presence of relic gravitons may have perturbed 
the early universe in ways that might be observable in the fine details of the 
CMB background. These gravitons would introduce a small dispersion in GWs, 
which might then leave an imprint on the CMB. We will discuss the potential 
presence of such an imprint in next Section. 

Now, let us expand the field 0(j according to harmonic oscillator operators 
b, , as a simple model of a string. The fields are expanded as 
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(j)i^^ = { — )Y^ E^b{k)e'^'^'''> + h^e""^^^ (21) 

where is a tetrad, which is discussed more below. The summation runs from 
{— oo, oo}. The product (p'^'p^a = <5/i^y is a harmonic oscillator operator 

(22) 

Hi) Ekk' El,bik)b{k')e'<^(-^^+'''(-^"> + 6t(fc)fot(fc')e-»e(fe)-*e(fc')^ 

The sum gives a delta function on k and fc' and the first term is the Hamiltonian, 
which after the use of a commutator the RHS term is 

r.(k.^.= C2)EkEi.bHk)h{k)+ (i)Efci^^.b(k)b(-k)+&t(k)&t(_k), (23) 

where the ZPE term has been dropped. The first RHS term is a familiar Hamil- 
tonian type of term, while the second term is similar to a squeeze operator in 
quantum optics |23| . 

The tetrad E"^' is the amplitude of the field. This plays a role similar to 
the minimal electric field E ~ y^hoj/Veo in box normalization |24| . A plausible 
choice for tetrad is then E'^ = a'uiS^, where a' is the string parameter and w 
the frequency. For a' ^ l/u this is a small term. 

The curvature in quantum modes is then 

(i)^ (A;.A;,(ii;2)f^i?2^-fcpfc,i?2^(£;2)^)(6t(fc)6(fc)+fe(fc)6(-fc)+fot(fc)^,t(_fc)), 

k 

(24) 

which is 0{a) in the scale parameter. In fluctuations of the curvature the metric 
is g ~ 6L/L, for 5L > Lp. The connection terms are of order F ~ SL/L^ and 
curvatures are R ^ 6L/L^. The wave vectors are k ^ 1/L and the scaling 
parameter is a'uj ~ SL/L. From a dimensional and scaling perspective this 
answer appears at least proximal. 

For the sake of simplicity, let us write the curvature tensor as 

R^,,p ~ (i) ^H^,<,p(fc)(6t(fc)6(fc) + b{k)b{-k) + b\k)b\-k)). (25) 

A: 

The second order term is formed from the total contraction on the Riemann 
tensor RaSfiuRap, it is 
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(b'<{k)bikp{k')b{k') + b\k)b{k){l){k')b{-k' + b\k')b\-k')) + {b{k')b{^k') + 
bHk'p{-k')p{k)b{k) + {b{k)b{-k)+b\k)b'^{-k)){b{k')b{-k')+W{k')W{-k'))y 

(26) 

This term is to 0{a'^) and contributes a term 0{a'^) to the Lagrangian. 

Consider the operator matrix operation R^i,crpR^'^'^'^\'ni) ■ The first term has 
the operator matrix elements 

6t(fc)&(fc)&t(fc')6(fc')|m) = b\k) Y,n \n){n\b(k)b'^ {k')b{k')\m) 

(27) 

= m{k')n{k)5„inSkk' 
where E„ |'^)('^| is a completeness sum and the momentum values assumed in 
the states |m) and \n). This contributes an energy-squared. A similar analysis 
for {m\R^,„pR^'"'P gives 

(m|6t(fc)6(fc)(6(fc')K-fc') + b^'ik'pi-k')) = m{k){b{k')b{-k') + b'^ {k')b\-k')) 

(28) 

and for R^^^pR^""'P\m) 

{b{k)b{-k) + b'^{k)b'^{-k))b\k')b{k')\m) = m{b{k)b{^k) + b^k)b\-k))\m). 

(29) 

The operators b{k)b{—k) + b\k)b^ {—k) form the squeeze operator [53] 

S = exp((i)(z*&(fc)-z&t(fc))), (30) 

where z* = z = i {{1 / 4)11 pi,crp{k)IlP'^ "'''). Hence the action phase due to the 
action e^^ contains a squeeze operator. The final operator term is more compli- 
cated. The operator terms b{k)b{—k)b{k)b{—k) and W {k)h^ {~k)b{k)b{~k) are 
evaluated by commuting operators and this leads to the square of number opera- 
tors n{k)n{~k). The terms b{k)b{~k)b{k)b{^k) and b\k)b\-k)b'^ {k)b\~k) are 
then a product of terms which represent a squeezed state operator. The squeeze 
operator S{z) acts upon the displacement operator D{a) — cxp(a6^ — a*b) so 
that S{z)D{a) ^ D{a)S{z), 

S{z)D{a) = exp[{z*b^ - z{b'^)^)/2] exp(afet _ a*b) = 

(31) 

exp[(z*62 _ z(6t)2)/2 + Q,fot _ a*b) exp[-(i)(z*a6t - za*b)], 
which effectively creates a modified displacement operator 

exp[{z*b^-z{b'')^)/2+ab^ -a*b)S{z)D{a) = cxp[-{-){z*ab'' -za*b)] = D{z*a). 

^ (32) 
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The action of the squeezed state operator on b is SbS^ = 6cosh(|z|) + 
sinh(|2;|), which is a BogoUubov transformed operator [26]. For a set of bosons, 
here hnear gravitons, with the same state there is then X^n "-"^"/v^l") states 
with the operator acting on this X^n '^"(°^)"/v^ acting on the vacuum. This 
operator may be formed from the S{z)D{a)S^ (z) for a smah and \z\ ^ \a\ with 

S{z)D{a)S^{z) ~ exp[-(i)(z*a6^ - za*b)], (33) 

and where we may then define z*a/4: — >■ a, and the Bogohubov transformation 
of the operator b—b^ constructs a displacement operator ^7\. In this way the 
term in the action describes the squeezed state operator which acts on the field 
raising and lowering operators to define a displacement operator for coherent 
states, which in the case of photons are laser states of light. 

We then evaluate a Wilson loop [15] W{4>'i) = exp(J' i(j)^^e'^dx''). In the 
path integral 

Z[cj),W] f V[(j}]We'^^'l'l (34) 



The infinitesimal shift in the field ^ + adjusts W] Z[cj) + 5(i), W] 
(W^) and the expansion is 



{W) ^ J D[(f>]Wi(f> + ,50)e*-5['>+*>l 

= {W)+J D[m{^ + ^e^^"')' 
where the invariance of the expectation gives 



(35) 



SW ^ ^ W6S _ ^ ^ SlnjW) ^ .5S _ ^ ^^^^ 
S(f> 5(f> 5(f) 6(j) 

This formula is only well defined for a polynomial function. So we make the 
following approximation. The loop is considered to be very small and in that 
way we can approximate the Wilson loop with 

1¥(0) = l + i0^,e^Ja;", (37) 
so that the functional derivative of W{(t)) is 

^^^e.j;6ix-x') (38) 
for a unit area. The solution is then 

Now, let us consider the second order expansion 
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which gives the result 



where by continuing the series this leads to 

{W-^) ^ ze'<*-->^"*"''(5(a; - x') = i{e'^^''^"^'=")S{x - x'). (42) 

The input of an expansion of the field (j) results in the expectation of an operator 
with the form of the displacement operator. 

It is now important to understand the form the fields in the expansion in 
D{a). The Wilson loop is a form of the Stokes' law [29] and 

- iln{W) = ^a^^^^e'^de"". (43) 

In vacuum the canonical h^^ and /ixx polarizations obey = D/ixx — 

[T7] . The longitudinal modes due to terms obeys [121 US] 

ahc = m'^hc, (44) 

where the mass is a topologically induced mass. The longitudinal he — (j>^ then 
defines the equation [15l [lU [20] 

0(f)^^ = m'^cj)^^ (45) 

where a Lorenz gauge sets terms with D<j> = [151 [SI HO]. This term plays a 
role similar to the Helmholtz potential in electromagnetism 

<J> = 7^/"rf^^#^> (46) 
^T^eo Jv |r-Fr 

but in the case of f{R) theories it results an effective potential through the 
identifications [151 HQI 



$^/'(i?) and ^^M^l^Rm (47) 
which give a Klein - Gordon equation for the effective $ scalar field [TS] [20] 

»=^. (48) 

The 4)fj,i, which physically contributes to the Wilson integral has a source term 
which is the topological mass. 
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3 Potential imprint in Cosmic Microwave Back- 
ground 

Wc recall that CMB is thermal radiation filling the observable Universe almost 
uniformly [3TJ [52]. Precise measurements of CMB are fundamental for cos- 
mology, because any viable proposed model of the Universe must explain this 
radiation. The CMB has a thermal black body spectrum at a temperature of 
~ 2.7 K. The peak of the spectrum results in the microwave range frequency 
of 160.2 GHz [211 Eg. 

At the present time, the best available data on CMB arise from the WMAP 
Mission |2H I22j. which used symmetric, rapid- multi- modulated scanning, rapid 
switching radiometers to minimize non-sky signal noise. WMAP works within 
the Solar System and to take into account weak potential effects on CMB by relic 
massive GWs we can use the weak field approximation (the linearized theory). 
In the linearized theory, the coordinate system in which the space-time is locally 
flat has to be used and the distance between any two points, say A and B, is 
given simply by the difference in their coordinates in the sense of Newtonian 
physics [17] . This frame is the proper reference frame of a local observer, which 
we assume to be located in A within the Solar System. In other words, we 
assume that the space-time within the Solar System is locally flat with respect 
to the global distribution of CMB. Our goal is to understand how relic massive 
GWs perturb the trajectories of CMB photons between A and B. The global 
effect results a particular gravitational lensing [30] due to relic massive GWs. 
Some clarifications are needed concerning this issue. In our linearized approach, 
gravitational lensing can be described in the local Lorentz frame perturbed by 
the first order post-Newtonian potential. Hence, one can define a refractive 
index 001 [31] 

n=H-2|y|. (49) 

In the usual Geometrical Optics^ the condition n > 1 implies that the light in 
a medium is slower than in vacuum |32j . Then, the effective speed of light in a 
gravitational field is expressed by (SO] [31] [32] 

w = -«l-2|y|. (50) 

n 

Thus, one can obtain the Shapiro delay [S^ by integrating over the optical path 
between the source and the observer: 

nohserver 

/ 2\V\dl. (51) 

J source 

The situation is analogous to the prism [32) . 
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3.1 Gravitational lensing in the direction of the propagat- 
ing gravity wave 

For the sake of simphcity, assume that A and B are both located in the direction 
of the propagating massive GW which we assume to be the z direction. By 
using the proper reference frame of a local observer the time coordinate xq is 
the proper time of the observer A and the spatial axes are centered in A. In the 
special case of zero acceleration and zero rotation the spatial coordinates Xj are 
the proper distances along the axes and the frame of the local observer reduces 
to a local Lorentz frame [T7]. The line element is [T7] 

ds^ = -{dx^f + S.jdx'dx^ + Oi\x^'^)dx°'dx'^. (52) 
The connection between Newtonian theory and linearized gravity is well known 

m 

300 = 1 + 21/, (53) 

where V is the Newtonian potential. Let us consider the interval for photons 
propagating along the z -axis 

ds'^ = goodt'^ + dz'^. (54) 

The condition for a null trajectory {ds ~ 0) gives the coordinate velocity of the 
photons 

d7 

vl^{-f = l + 2V{t,z), (55) 
which to first order is well approximated by 

Vp^[l + Vit,z)]. (56) 

Knowing the coordinate velocity of the photon, the propagation time for its 
traveling between A and B, which corresponds to the proper distance AB in 
presence of the graviton, can be defined: 

Ti{t) ^ r -^T- r V{t', z)dz, (57) 

JzA -'0 

where T represents the uniform propagation time of the photon between A 
and B (i.e the proper distance between A and B in natural units) as if it were 
moving in a flat space-time, i.e. in absence of GW, and t' is the delay time 
which corresponds to the unperturbed photon trajectory: 

t' =t-{T-z) (58) 

(i.e. t is the time at which the photon arrives in the position T, 
soT - z = t-t'). 
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In order to compute Ti we need to know the Newtonian potential V{t, z) which 
is generated by the massive GW. We recaU that the effect of the gravitational 
force on test masses is described by the equation 

X' = -Ri,ko^\ (59) 

which is the equation for geodesic deviation in this frame [T7], and i?ofco 
linearized Riemann tensor jl7) . 

On the other hand, with an opportune choice of the Lorentz gauge, the 
linearization process of f{R) theories which generates the third longitudinal 
mode he — hc{t — vgz) enables a conformaUy flat hue element [TSl [TBI [20] 

ds^ = [1 + Kit - VGz)]{-dt'^ + dz^ + dx^ + rfy^), (60) 

where vg represents the group velocity of the massive GW. In fact, the velocity 
of every standard massless tcnsorial mode /i^^ is the light speed c, but the 
dispersion law for the modes of he is that of a massive field which is a wave- 
packet [HI [ini HOI ■ Also, the group- velocity of a wave-packet of he centered in 
^ is [15l[T6l[20] 

(61) 

which is exactly the velocity of a massive particle with mass m (see Eq. (jH])) 
and momentum jt. 

This group-velocity is function of both of the mass and frequency of the 
wave-packet [HI HH ED] 



VG = ■ (62) 

Even if the coordinates are different from the coordinates we recall 
that the linearized Riemann tensor is gauge invariant [17] . Hence, we can 
calculate it directly from Eq. ((501) • Following [T^j it is: 

RfiiyalS = -^{df^dpha^ + d^dah^p — dadfjh^^ — d^d^ha/j}, (63) 
that, in the case eq. begins [IB] 

-^070 = ^{d^doherjoj + dod^heS^ - d'^djhe'noo - dodoheS"}; (64) 
the different elements are (only the non zero ones will be written) [16] : 

dfhe for a = J = 
d^dohevoy = { } (65) 

-dzdthe for Q! = 3;7 = 

d^he for a = 7 = 
dQd^he6^={ ) (66) 

dtdzhe for Q! = 0;7 = 3 
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—dfhc for a ~ J ~ 
d^hc for a = 7 = 3 
-dtdzhc for a = 0; 7 = 3 
dzdthc for a = 3; 7 = 
d^hc for a = 7 . 



By putting these results in Eq. ([M]) one gets 



-^010 ~ 



d2 — _lh 
^010 ~ 2 



^030 



Let us put Eq. dH]) in the third of Eqs. We obtain [H] 

1 



030 — 2^ 



which shows that the field is not transversal. 
In fact, Eq. implies [IS] 

2 



and 



-rn^hcit — vgz)z. 



(67) 



(68) 



(69) 



(70) 

(71) 
(72) 

(73) 



Therefore the effect of the mass is exactly the generation of a longitudinal force 
(in addition to the transverse one). Note that in the limit 7ti — > the longitudinal 
force vanishes. 

Equivalently we can say that there is a gravitational potential |161 117) : 



1 1 

V{l^,t)^--h^{x^ +y^) + -m^ J h^{t ^ VGa)ada, (74) 

which generates the tidal forces, and that the motion of the test mass is governed 
by the Newtonian equation [1^1 [T7] 



(75) 
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Now, we can use Eq. dTJ) to compute Ti in Eq. dST]). Since = V{x,t) + V{y,t) + V{z,t) 

we get 



Ti{t)^T- [ Viz, t')dz = T - l-m^ [ [ ^c{t' - vaojada. 
Jo 2 Jo Jo 



(76) 



Thus, the variation of the proper distance between A and B from its unper- 
turbed value T which is due by the presence of the massive GW he is 

STi{t) sa im^ dz hc{t — T + a — vca)ada = 

= jm^ hc{t — vqz — T + z)dz — jm^ h'^[t — T + a — VGa)z^dadz. 

(77) 

Introducing the Fourier transform of he defined by 

/oo 
dthc{t)cxp{iu!t), (78) 

eq. (|77p can be integrated in the frequency domain by using the Fourier 
translation and derivation theorems 

^ = (79) 

where 

= i^'i^^){^-p^^nvG - 1) - 1+ 

+ [T^cxpiioTivG - 1) - 2Texpiwr(wG - 1) + 2cxpiujT{vG - 1) - l] - 

(80) 

is the longitudinal response function for relic gravitons. 

In order to use eqs. (|79p and ((80)) we recall that relic gravitons represent a 
stochastic background [HI [20]. Hence, one has to use average quantities [T8l[20] . 
The well known equation for the characteristic amplitude [TH], adapted for the 
third component of GWs can be used [50] : 

h,,{f) ~ 1.26 X io-i«(^)y';i2^o^^9-(/), (81) 

obtaining, for example at 100 HZ and taking into account the bound ([20]) . 



h^ciWOHz) ~ 1.7 X 10"^^. (82) 

Considering a graviton propagating with a speed of vq = 0.999 (ultra- 
relativistic case), if we insert these values in eqs. (j79p and (|80)) we get T(w) « 
0.02 and 5Ti « 3.4 x lO^^^m for a proper distance between A and B of unper- 
turbed value T = 1km. The situation is different for a speed of 0.9 (relativistic 
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case). In that case one has T{uj) sa 0.19 and 5Ti 3.4 x 10~^^to. For a speed 
of 0.1c (non relativistic case) we have T(cj) sa 0.99 and 6Ti sa 1.6 x 10~^^m. 
The situation is better at lower frequencies. For / = IQHz eq. ([5T|) gives 
hcc — 1-7 X 10~^^. The response functions result practically unchanged, there- 
fore wc gain an order of magnitude, i.e. STi « 3.4 x 10~^^to for vq = 0.999, 
Sfi sa 3.4 X 10-23to for vq = 0.9, and Sfi sa 1.6 x W'^'^m for vq = 0.1. 

Here wc discussed the variation of the photons' paths in the z direction which 
is the direction of the propagating relic GW. Clearly, analogous effects, which are 
due by the transverse effect of the GW (eqs. ((7T|) and ([7^ ). are present in the 
X and y directions. Thus, eqs. ([71]) and ([50)) can be used to discuss the general 
gravitational lensing in our model. We developed the complete computation in 
the z direction, the extension to the x and y directions is similar. 

The global effect of these variations of the photons' paths in CMB should be 
analogous to the effect of water waves, which, in focusing light, create optical 
caustics which are commonly seen on the bottom of swimming pools. 

4 Chaos and Relativity in Orbital and Optical 
Systems 

The consequences of GWs form /(i?) theories arc observable fingerprints on 
the structure of the universe. Massive GWs will act as lenses which generate 
caustics in the motion of light and other particle fields. These caustics will then 
have measurable influences on the CMB or upon the distribution of galaxies in 
the universe out to z = 1 and beyond. The following looks at the issue of how 
general relativity can amplify chaotic dynamics, and further can amplify optical 
chaos. This is illustrated in a three body problem, and in an elementary optical 
model. This digression into another aspect of relativity is meant as a way to set 
up analysis for the phenomenology of massive gravity waves. This illustrates 
how to proceed through the examination of elementary systems. The extension 
to more complex structures, such as a many body problem of galaxies and dark 
matter, will require numerical methods. This indicates that relic gravitons could 
have a larger effect on the early distribution of matter. 

One of the early tests of general relativity was the prediction of the perihe- 
lion precession in the orbit of Mercury |17j . This is a departure from Newtonian 
gravity that is largely post- Newtonian to 0(l/c^). These general relativistic cor- 
rections are completely integrable and with no chaotic dynamics associated with 
them. In a three body problem, with a large central mass, a larger distant mass 
which is treated as Newtonian and a smaller satellite with 0(l/c^) relativistic 
departures, will exhibit chaotic dynamics in the small body. The additional rel- 
ativistic corrections will interplay with the irregular chaotic dynamics and are 
shown below to contribute to a Lyapunov exponent [34]. In effect a Lyapunov 
exponent A = log{K) will have a relativistic correction A = Aq -I- A(0(c~^)), 
and this correction then amplifies the chaotic behavior of the system. This is 
extended to optical systems. Einstein lenses [35] are a Newtonian gravitational 
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phenomenon, and general relativistic corrections to 0(\/(?) are minor, for the 
impact parameter on such a gravitating body is too small to be observationally 
significant. Yet for a complex Einstein lens, say analogous to a compound lens 
due to smaller scale clumping of matter, a light ray may have a succession of 
small angular deviations. These angles of deviation will have a compounding 
effect similar to angle deviations of a particle in an arena. This will result in 
increasingly complex optical caustics which in analogue with chaos are difficult 
to predict. This is further compounded if the gravitating clumps of matter are 
difficult to observe directly, such as with dark matter |36) . In a manner sim- 
ilar to the case with orbital dynamics general relativistic corrections may also 
enhance this optical chaos or turbulence. This Section connects two different 
aspects of chaos and relativity to present issues with the analysis of three body 
systems with parameterized post-Newtonian parameters. Subtle enhancements 
of chaotic dynamics or the increase in a Lyapunov exponent might be docu- 
mented in such a system. This should then be an observable characteristic of 
complex relativistic systems. The optical analogue illustrates how fine detailed 
structure in a distribution of matter which is an Einstein lens could influence 
the complexity of optical caustics. Localized regions of large gravity fields could 
then further amplify this complexity as well. This might lead to methods for 
mapping any local density variation in dark matter. It may further lead to an 
understanding of how a stochastic distribution of massive gravitons can lead to 
the distribution of galaxies in filaments and walls. 



4.1 General relativity to 0(l/c^) 

In general relativity the equation of motion for a test mass particle around a 
fixed central mass is [17] 

(fu GM SGMu"^ 

Here I is the constant specific angular momentum. We recognize this dif- 
ferential as the harmonic oscillator equation of Newtonian mechanics with a 
constant force GM/P, plus the term ^ {u/cf'. The anomaly angle obeys the 
dynamical equation |17j 



dO _ I 
ds 

and 



lu"^ (84) 



dt _ E 
Ts ^ 1 - 2GMuld^ ^ ' 

for E the potential energy per unit mass of the particle " at infinity" = constant. 
For GMjf? ^ 1 we may solve this problem by perturbation methods. The 
solution of interest is 0(1) plus 0(c~^), which would be Newton plus first order 
GR correction. The expansion is carried out with the variables u, according 
to 
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(86) 



e = eo + eei+0{e^). (87) 

Here the term e = 1/c^ gives the order of the expansion. The differential with 
respect to 9 to first order in e is taken as 

d d d . 

d0 ~ de^ ^ '^d9^' ^ ^ 



If we input the expansion for u in equation ()86p into the differential equation of 
motion ()83p the following two equations are obtained: 



: ^ + "1 - Skuq = 0. (90) 
d^o 

The term d?uo/d9od9i = since uq is not a function of 6'i. Further, the term 
K = GM . The 0(1) differential equation ((5^ has the solution 

uo = ^(l + e'co.s(6'o+a)), (91) 

which is the standard Newtonian solution for the radial velocity for a particle 
with orbital eccentricity e' and anomaly angle a |17] . Now, let us consider on 
the expansion of 9. We set = 1 and insert this into the equation for the 
angular velocity equation 

d9 

= {'^- 2Keu)lu'^, (92) 

where 1 — 2Kew is the Schwarzschild transformation between proper and standard 
time coordinates |17| . This differential equation has the two contributing parts: 

0(1) --^-lul (93) 

0{e) : = 2luo - 2kIuI (94) 

We are primarily concerned at this point in the solution to order 0{e) for the 
orbit of a test mass in a GR orbit, 



d^m 



d9, 



ui - 3km5 = (95) 



where the Newtonian solution uq is given by equation (1911) . The square of uq in 
the non-homogenous term is 

ul = (^^y(l + 2e'cos{0o + a) + e'^cos^{9o + a)y (96) 
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which by elementary trigonometric identities is 

ul = (^) ' ((1 - e') + 2eW((0o + «)/2) + e'2cos2(0 + a)) . (97) 

The solution is elementary at this point. The first non- homogeneous term is 
going to give a solution 

~ ■^(l-e')(l + e'cos(0o + «)), (98) 
and the quadratic trigonometric functions determine the solution: 

Ul = ^(^{l+e'cosieo+a)) + ^-{cosi9o+a)-3) + '-^{cosi2ieo+a))-3)y (99) 

The hard part is the perturbation of the third planet. The Jovian planet obeys 
a similar dynamical equation, but where c — >■ oo and Newtonian dynamics is 
recovered as 

Here v = l/r2 for this additional planet, and we define u = uq + eui ~ l/ri. 
Similarly, the angular momentum is defined by |17j 

do' _L o , ^ 

The angle 9' may exist in a different plane than 6, yet as an approximation we 
put both angles in the same plane of motion. Now we need the coupling between 
the two bodies. We assume they are Newtonian as 

F GMm ^l^ "I^,, , (102) 
|ri - r2|'^ 

which is approximately 

GMm I ^ , 3 r 1 • f 2 

'2 



To find the distance jri — ^2! we consider the plane of the two orbits as complex 
valued and that the positions of the test mass and the larger mass are give by 
f\ = rie*^i and r2 = r2e'-^'^ and so the distance between the two masses is given 

by 

- r2|' ^rl+rl~ 2r2r2cos(0i + 62). (104) 

The potential energy 

Uin,r2)^-^^ (105) 
|ri - r2| 

defines the force in equation (|102p by F = — VC/. For r2 ^ ri the denominator 
in the potential may then be cast in the u, v variables 

U{ui,U2) ~GMmv(l- (^^y + 2^cos{{uji + uj2)t)y (106) 
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Here uJi = d9i/dt, for i = 1, 2 for the two bodies. This is the perturbing potential 
for the two orbits of the bodies in the same plane. 
The total Hamiltonian is then 



H = +eH 



+ H'i° + n'v[ (1 - (^) ' + (^) cos{e + 9')] 



(107) 



where the first three terms are harmonic oscillator Hamiltonians 

Uh° — 1 ,i2 I 1 ^ ,,2 TTho _ 1 2 I 1 ,2 TTho _ ^ ^2 , 1 '^'^ ,2 n nS^ 

-^0 - 2^^o + 2^"0' ^1 - 2^*1 + "i' -2^+2L2 ^ ' 

We now have two order parameters e = 1/c^, and another k' = Gm, where the 
mass m is the mass of the "Jovian" planet. The Hamiltonian term that scales 
according to ek' for 6* = 6*0 + eOi is 

^ -.'(l - ^) (^)\r(l - ^) (109) 

To compute the Lypunov exponent explicitly the gradients of the Hamiltonian 
with Pij Pv and uq, ui, v are first found. With v/uq <C 1, Ui ^ uq these 
are then to order [v/u^Y 

VpgH^iio, Vp^H = eui, Vp^H = Uy (110) 



V„„i? = t^tUq ~ GeuQUi — k' ^cosiO + 6*') (HI) 
/2 u;^ 



Vu^H ^ —ui ~ 3euQ ~ eK — (112) 



V^H = ^v + k'(i-2^ + 2—cos{e + 9')]. (113) 

These are the forces F ~ ~VH due to the three configuration variables mq; ui 
and u. The last right hand side terms in V„ji7 are dependent upon both the 
general relativistic correction 0(l/c^) and the gravitational coupling with the 
Jovian planet k! . 

We consider the change in the phase space flow 

Z + = (m + Au, p + Ap) (114) 

The change in momenta due to the perturbation from the Jovian planet is 

2 2 2 

Ap ~ m\k'(i -2^ + 2—cos{e + B')] - K'^cos{e + 0') - ek'^I , (115) 
I \ ui uo y ui uil 
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where the last term is a coupling of general relativistic 0(l/c^) effects and 
planetary perturbation. To 0{k' /c^) Au oc Ap. Define Ap{t) to be the deviation 
in momentum due to planetary perturbation, and let 6p{t) be the deviation due 
to the 0(l/c^) coupling term. The Lyapunov exponent is then 

A . hm i,4M^)±M) ) ^ 1,, i U^EK) + MMM^l (116) 



so that 



Ap(to) I *->oo i L \Ap{to)J Ap{t) 

A-A„ + .(^)\ (117) 

with Ao defined for e = 0. The exponential divergence in phase space be- 
tween nearby trajectories has then contribution in addition to Aq with Zit) ~ 
Z(tQ)e'^°*e^^'"^^°^ *. Thus general relativity will bring about the onset of chaotic 
behavior, or the breakdown of numerical unpredictability, earlier. 

This should then manifest itself in semi-relativistic systems with three bod- 
ies. A system, such as two neutron stars in a mutual relativistic orbit with 
a third companion further away and executing Newtonian dynamics, of this 
sort will then have more chaotic behavior which is amplified by general rel- 
ativistic effects. This simplified model suggests that a general parameterized 
post-Newtonian-multibody perturbative theory is needed. Such a model will 
then be more suited for the examination of complex general relativistic systems 
that include several bodies. 

4.2 0(l/c^) Optical Corrections in Einstein Lensing 

The Einstein lensing of light is now a common observational feature of deep 
space astronomy since the launch and repair of the Hubble Space Telescope, see 
[37] and references within. Here a complex optical gravitational lensing system 
is discussed with some analogues to the mechanics above. A large elliptical 
galaxy will have an overall gravitational lensing effect, but there may be sub- 
lensing as well if the density of dark matter exists has some variation. This 
results in a type of optical turbulence, analogous to chaos. Further, this may 
also be amplified by general relativistic effects. Unknown configurations might 
exist with dark matter density increasing in the vicinity of a large black hole. 

The general theory of gravitational lensing [30l [31] shows that a light ray 
which approaches within a radius r ^ 2GM/(? will be deflected approximately 
by an angle 6 = GM/rc?. In a more general setting the deflection of light is 
given by the Einstein angular radius 



where dis, di, ds are the angular diameters to the gravitational lens, the source 
and the distance between the gravitational lens and the source. For dis, di, ds 
the angular diameters to the gravitational lens, the source and the distance 
between the gravitational lens and the source. The condition da = di + dsi 
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obtains locally where cosmological frame dragging is small. This theory is the 
weak gravitational lensing approximation, where the deflection of light is essen- 
tially a Newtonian result [30] ■ The distance relationships are determined by 
Ods = fids + a'dis- The reduced angle of deflection a{9) — {dis/ds)a'{9) gives a 
relationship between the angles of importance a{9) + (5 ~ 9. 

Complex distributions of matter can act similar to a compound lens in a weak 
gravitational limit. However, light rays which pass close to clumps of matter to 
exhibit 0(l/c^) deviations will exhibit deviations from this linear summation. 
A light ray which passes through a set of random lenses will display caustics 
which are similar to the caustics seen on the bottom of a swimming pool. Fine 
grained structure in an Einstein lens can exhibit caustics which occur due to 
nonlinear perturbation in the density profile of matter. This nonlinearity in the 
symmetry of the lens will produce caustics which are analogous to chaos. The 
occurrence of a caustic has its connections with catastrophe theory |38] and the 
onset of a fold, which is also a mechanism for the bifurcation of vector fields in 
Hamiltonian chaos. 

For the position of a source x, the propagation of light along the z axis from 
this source then reduces the visual appearance of the object to ^ = £.y) 
along the axis of optical propagation. The weak gravitational lensing of light 
[50] then indicates that the deflection of the appearance of this object along the 
axis of optical propagation is given by 

Ae = V$(0, (119) 

for ^ the position of the image with the mass present and '&(^) the gravitational 
potential. The difference in the vector position of the image — is the 
difference between the position with the mass present and without it being 
present. The potential term obeys the Poisson equation [T3] so that 

V^$ = 2— ^ (120) 

The integration over the direction of propagation then gives the mass density 
in the plane of the image, often called the surface mass density S(^). The angle 
of deflection a is then determined by the Poisson equation and the potential as 

4G f{^-CMC),2,, 



for S(^) a mass/area density distribution in the image. The function S(^) plays 
the role of an index of refraction based upon a mass distribution, which for 
a thin lens will give the angle of deviation. For a gravitational thin lens, a 
weak field that is very small compared to the optical path length, and S(^) is a 
constant. The deflection angle is simply 

a{0 = ^3-^ (122) 



23 



where for small angles \^\ = £, = di9 and 



for the critical mass density Sc = {c^ /ATTG){ds/disdi). This is the minimal 
mass density which might be distributed in the area of an Einstein ring [39] , 
For a more complex arrangements of gravitational lenses, such as large density 
nonlinearities, the mass density S has a general form 

= JdzaidiC:,, diCu z) (124) 

The position of the image then plays the role of the vector r in equations 
and beyond in the above discussion. Let the reciprocal of the vector 
norm |^| = Xjv plays the role of u. The analogue of the Newtonian equation of 
motion in equation 1 with c — > oo is then 

^^=J%^ (125) 

for ^ the impact parameter. Now the Newtonian description of gravitational lens 
deflection has the effective photon angular momentum per mass term j = (c^)~^. 
The general relativistic extension of this equation is then 

To order expansions the analogue of equations (|89| and (|90|) are 

a™ = (127) 

A.. = (128) 

for uq the reciprocal of ^. The term k = GM for a general distribution in the 
plane of the Einstein ring is k/^ ~ AirGYidisdi/ds, which reproduces the Einstein 
ring case in the first approximation. The second order term is Av = A^/^^ and 
a ~ 2A^/^, which reproduces the weak field gravity lens result. The 0(l/c^) 
correction gives an effective general relativistic correction term 

Aa.-i-(-) (129) 

This correction term is not likely to be detected directly by extra-galactic 
sources, such as the dark matter lensing of light by the Abell galaxy cluster 

Bol. 



24 



4.3 Optical Chaos 

Just as the 0(l/c^) correction to Newtonian dynamics enhanced chaotic dy- 
namics, or contributed to a Lyapunov exponent, we might expect a similar 
amplification of optical chaos, or the statistical appearance of caustics by the 
clumping of matter. Small local region where gravitating mass is clumped to- 
gether will result in the deviation of the light ray by some small angle 60, which 
is an error in computing the subsequent tracing of the ray. With a succession 
of n such small deviations the first angle deviation is amplified by ~ 2"i56'i, the 
next by ~ 2""''^(56'2, where for large n and Oi ~ Mi the total angular error in 
computing a ray trace will be approximately 2""'"^(50. This is analogous to the 
arena problem of computing the trajectory of a ball. 

The vector ^ describes the visual appearance of a distant object along the 
axis of propagation. This vector describes the deformation of a wave front by 
the lensing action of the intervening gravitating body. The gravitational lens 
is usually considered as a symmetric lens [50]: but nature may provide local 
clumping of material which introduce some chaos in the ray tracing. Further, 
the over all gravitating lens may be sufficient enough to produce small 0(l/c^) 
relativistic deviations from a purely Newtonian lensing. Above the formula 
for this relativistic deviation is given. What is then needed is an analogue to 
Lyapunov exponent for the classical unpredictability of a ray trace due to New- 
tonian gravitational sources. A multiple set of ray tracings is then a description 
of the deformation of an electromagnetic wave front, and perturbations on the 
vector ^. In what follows such a development is presented to describe the chaotic 
perturbation of this vector. 

The propagation of a plane electromagnetic wave front is given by tp{r) = 
■i/'oe''^ ''"'^* . The occurrence of a gravitational lens perturbs the the wave front 
according to 

^'{r) = (x(f)e^^«)^.(f). (130) 

Here the is the change in the wave front phase and x(r) is the change in the 
wave front amplitude. The vectors describing the visual appearance of the image 
are ^ = Vr^4'{'r) for r\\ coordinate directions along the wave front. This means 
that A(/)(^) = $(^), which is a Poisson equation [13]. The phase deviations are 
caused by an effective index of refraction [301 131] in the Newtonian limit, and 
the gravitational potential is the source in the Poisson equation. A Gaussian 
random distribution of sources results in the second order structure function 

^0(/?) = (l0(r)-0(r-p)P). (131) 

The vector p = ^ + where ^ = r — /?, so D^{p) is a phase variance between 
two different direction in the aperture plane. 

The phase terms obey a Poisson equation, where some distribution of sources 
is present. For optical perturbations compatible with the second order structure 
function the gravitating perturbations arc in a Gaussian distribution, where 
Gaussian distributions of perturbing sources means that the equation ()120p 
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becomes 

y2^^2^ l^^TT «?/2a^ (132) 

47r (y2^a)3 ^ ' 
where each <C ■ Each one of these sources gives a solution 



= {^i|A^,){l|r)Erf{r/^/2cJ), (133) 

for the variable = r, and the solution converges to a point course in the 
limit CT — >■ 0. Each of these perturbing changes on the aperture vector is due to 
a succession of matter clumps. A photon which passes close to each clump is 
modeled as having its angle deviated, and its path is then stochastically deviated 
away from a path given by equation (jl32p . The small angle of deviation for 
a(^) — > d(^) + (5d(^) is determined by the Gaussian distribution as 

for p{^i) = (l/47r)(/^/(-\/27rcr)'^)e~^' /'■^'^ For simplicity the angle of deviation 
(5(3(^) will now be treated as a scalar and with ^ the angle deviation is 

5a^K f\[C^p{S,[)dS,[ (135) 
•' i 

such that a ~ (C^^)- This is a partition function analogous to that in the 
Ising model [H], but here instead of a set of spins that exist in space there are 
stochastic angle changes in a ray trace of light. In this particular model these 
stochastic angle changes are assumed to be on average the same. 

This partition function can be demonstrated to be similar to the Ising model. 
For the variation in the stochastic variable = — Cj-i iii the exponent, the 
product of any two variations vanish SS^jS^j ~ 0, so that 

for i = j the sum of these stochastic variables is 

, n-l (n-l)/2 
3=0 j=l 

This means there exist additional endpoint terms which do not conform to the 
Ising type of construction. However, for a large enough n this error should be 
minimal. The expectation is approximately 



^/ n rfe,r'e.Tp(-e,_iO/?), (138) 
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for /3 = 1/tT^. /3 is analogous to the Boltzmann factor, which is determined by 
the scale at which matter is lumped together. A correlation length scale is 

~ l/log{tanh/3), (139) 

which for /3 ^ 1, or equivalently for large cr is A ~ cr. This approximate 
formula is a ray trace path analogue of the Lyapunov exponent in time, which 
determines a length A where the prediction of a ray trace breaks down. This also 
illustrates that this breakdown of ray tracing occurs on a scale comparable to 
the length scale of the perturbing. This loss of ray trace prediction is manifested 
deformations of the angle deviation across the aperture distance, or deviations 
in the symmetry of an Einstein ring. 

The goal now is to determine if there are enhancements of optical chaos, 
analogous to optical turbulence in the Earth's atmosphere, due to 0(l/c^) cor- 
rections. To examine this we consider the angle of deviation due to Newtonian 
gravity, optical path length turbulence, and relativistic corrections as conjugate 
to action variables J, J' , J" , and the path is given in a classical setting by the 
action 

{H + H' + H")dt = {JdaiO + J'da'{5£,) + J"da"{Ci). (140) 

For the angular momentum variables J ^ J' ^ J" in this equation then the 
action is entirely governed by the angle deviation, which for da = {da/dt)dt 
expresses this as a principle of least time. Just as in the case of planetary 
motion. The angle of deviation due to dumpiness of matter is approximated as 

Sac — Kexp{£_'^log{n)l3), (141) 

for n regions of matter or dark matter clumping. The region where the light 
ray is the most distorted by gravitating bodies is of a distance ~ ntr = ^ /?/2, 
which then gives an approximate relativistic 0(l/c^) correction 

^a,.5^^!^^!^^!±^(^)^ (142) 

where SO ~ Sac- An approximate Lyapunov exponent is then 



A^ lim izogfl + K(i^)[e«''°^'(")^+3«na2(02^20e«''°9("'^)f- 

n -¥ oo n V OTic L 

(143) 

Here there is an amplification of the ray trace uncertainty, or chaos, by the 
introduction of 0(l/c2) term as seen in the term 26'e«°'°^'(")^(S/i;c)2. For Ja^ ~ 
9 the contributions to the chaotic ray traced path from relativistic corrections 
and chaos are comparable and will contribute equally to the randomness of the 
caustic gravitational lens. 

The difference in the perturbed aperture vectors AS^ = S£^i — d£,s^s^^ de- 
termines the magnification M = d{S£,i)/d{S£^s)- From Hamilton's equations this 
is generalized to 

ASi= W^H ^ 2K^log{n)dfac{l + ^^{yY) (144) 
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with the deviation magnification computed accordingly. For this written ac- 
cording to the radius of curvature i? of a surface for a ray curve along a line of 
sight we have that 



The magnification for 77 = 2tt ^^2\R\/Yjc is = 1 + 778 + 0{'q^). The curvature 
R defines a tangent for the ray trace, which defines a caustic when the line of 
sight is along this tangent. The caustics along lines of site occur at swallowtail 
folds in the magnification map. 

For the orbits in different planes the above must be generalized some. Simi- 
larly the angular components arc given by the tangent vector parallel to pi and 
from there we may find the vector to each dark matter clump in its plane of 
motion with coordinates {u,6} and Further, as the angular momentum 

vector is given by Z = ri A pi/m (similar for L), then L is rotated relative to I 
by the Euler angles a, /3 and 7. The rotation matrix is then 



With these we may be able to put the problem in a general setting. This part 
is yet to be worked, and their may be resources to aid in this effort. 

The formation of filaments and domain walls of galaxies is then proposed 
to occur by this mechanism. The massive gravity waves in the very early uni- 
verse, such as in the post inflationary period, deviate the motion of relativistic 
particles in a manner similar to the optical focusing of light. These focal points 
of matter then set up their own gravity fields which persist through the sub- 
sequent expansion of the universe. A mesh of caustics with swallowtail cusps 
heuristically may be seen to produce a web of regions where mass-energy is con- 
centrated. The distribution of dark matter may then be established by caustics 
of gravitons and gravity waves in the early universe. 

5 Conclusion remarks 

In this paper the production of massive relic coherent gravitons in in a particular 
class of /(i?) gravity which arises from string theory and their possible imprint 
in CMB have been discussed. The key point is that in the very early universe 
these relic gravitons could have acted as slow gravity waves. They may have 
then acted to focus the geodesies of radiation and matter. Therefore, their 
imprint on the later evolution of the universe could appear as filaments and 
domain wall in the Universe today. In that case, the effect on CMB should be 
analogous to the effect of water waves, which, in focusing light, create optical 
caustics which are commonly seen on the bottom of swimming pools. This issue 
has been carefully analyzed by showing gravitational lensing by relic GWs, i.e. 
how relic massive GWs perturb the trajectories of CMB photons. 

The consequence of the type of physics discussed has been outlined from the 
point of view of an amplification of what could be called optical chaos. 




(145) 



[R] = [cos(7)] [1] [cos{a)\ . 



(146) 
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